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Abstract 

The effect of instability on the propagation of solitary waves along one-dimensional discrete non- 
linear Schrodinger equation with cubic nonlinearity is revisited. A self-contained quasicontinuum 
approximation is developed to derive closed-form expressions for small-amplitude solitary waves. 
The notion that the existence of nonlinear solitary waves in discrete systems is a signature of the 
modulation instability is used. With the help of this notion we conjecture that instability effects on 
moving solitons can be qualitative estimated from the analytical solutions. Results from numerical 
simulations are presented to support this conjecture. 
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I. INTRODUCTION 



The discrete nonlinear Schrodinger equation (DNLSE) is one of the most investigated 
systems in dynamical nonlinear lattices. In particular the DNLSE with cubic nonlinearity is 
an ubiquitous dynamical-lattice system which has been extensively studied because its direct 
physical applications, such as Bose-Einstein condensates (BECs) in deep optical lattices or 
optical beams in waveguide arrays, among others. In particular, the transport properties of 
neutra. atoms in EEC arrays has gained interest in the .ast few years Q due, in prineip.e, 
to the possible technological applications as matter-wave interferometry [2| or quantum 
information processing 3j, \4\. 

Of course, solitary waves are promising candidates for the coherent matter wave trans- 
port in BEC arrays. Notice that, in the strict sense, solitons (defined as solitary waves that 
maintain their shape and amplitude owing to a self-stabilization against diffraction through 
a nonlinear interaction) do not exist in the DNLSE. It is because the DNLSE is a nonin- 
tegrable system. However, approximate moving solitary-wave solutions can be analytically 
calculated. One of the main characteristic of these waves is that they emit radiation when 
moving, and eventually become pinned by the lattice This is a radiative deceleration 
which is extremely slow 5J], and therefore do not play a role in the present analysis. 

In the following we shall use a weak definition of soliton to refer to solitary wave solutions 
following from approximate integration of the DNLSE. 

We note that diverse analytical approximations have been developed in the last two 
decades to obtain moving solitary wave solutions of the DNLSE for bright solitons (the 



theory of the perturbed Ab 
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owitz-Ladik equation 



discrete approximation [121 ]. variational approximation 
theory 



construction 
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IOH . iteration methods 



ll| , semi- 



131 ] . or an explicit perturbation 
) as well as for dark solitons (quasicontinuum approximation 15j, or by 



17|). Of course, the DNLSE in the continuum approximation leads to the 



standard nonlinear Schrodinger equation (NLSE). Notice, however, that the soliton solutions 
of the NLSE do not move along the lattice, i.e. they behave as breathers in the DNLSE 
system. 

In the absence of dissipation and inhomogeneities 

n 

account that radiative deceleration is very slow [5j the main effect which reduces the life 
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duration of these solitary waves is the instability of the system. The instability effect in 
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moving solitons manifests itself as a distortion of exponential-like nature of the envelope 
during propagation |l4| . 

It is a common notion in discrete systems that solitons exist or can be excited in the region 

nnnn 

of parameters where modulation instability (MI) of planewaves appears [7|, l22j, |23j, l24j . So, 



solitons by themselves can be considered as a signature of the MI. Moreover, since the MI is 
a permanent feature of the discrete system, i.e. it does not disappear after the appearance 
of solitons, it is reasonable to conjecture that the continuous deformation of moving solitons 
and breathers after formation is a MI effect. 

The usual analytical procedure for stud ying M I effect on lattices is performing modulation 
stability analysis of planewaves 
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24, 



251 ] . From this analysis usually qualitative 



conclusions are derived for solitons. Recently the stability of narrow discrete modes in the 
DNLSE was studied for the case JU < in Eq. (HJ) |26|. 

Since solitons can be considered a signature of the MI effect in discrete systems (3, 0, 
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241 ] . it is natural to conjecture that approximate analytical forms of the soliton solutions 



may contain already qualitative information of the MI (e.g. strength and parameter region 
of existence). Here, the MI strength refers to the grade of how during the time evolution 
the MI effect distorts the initial analytical soliton solution. 

In the present study we are interested in studying anew the the stability of moving 
solitons in the DNLSE with the help of approximate analytical soliton solutions. In order to 
find these analytical soliton solutions we use a self-contained quasicontinuum approximation 
(SCQCA). Our aim is to show that these approximate analytical soliton solutions can be 
used to derive qualitative conclusions of the soliton stability in the range of parameters 
where the solutions are valid. In this regard we show that the amplitudes of these soliton 
solutions define an analytical upper boundary for a "self-defocusing" instability not only for 
bright but also for the dark solitons. Here the term "self-defocusing" instability refers to 
a decaying amplitude in moving solitons which also undergo width broadening. Numerical 
results supporting the analysis are presented. The standard modulation stability analysis is 
also revisited for comparison issues. 

We want to comment that it may sound superfluous to propose another method for 
calculating solitary waves in the DNLSE since other powerful methods as, e.g. those propose 
in Refs. 0,0, have already captured in fine detail interesting aspects of this equation. 
However, in contrast to those methods, the SCQCA can be straightforwardly extended to 
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two 31] and three spatial dimensions, where no other analytical methods for moving solitons 
have been proposed so far. Therefore, understanding the effect of the modulation instability 
on the solutions presented here is important in order to proceed to higher dimensions. 



II. THE SCQCA 



the DNLSE 



The DNLSE with cubic nonlinearity, which is the basic model for the one-dimensional 
BEC arrays in the tight binding limit pi 
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2l| . reads 



z3tV>n(*) + J{i>n-i{t) + Vwi(O) - u\^ n {t)\ 2 ^ n {t) = 0, 



(1) 



where ip n is a complex amplitude of t 



to the tunneling rate 
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re BEC mean field at the site n, J is proportional 



2l| and the nonlinear coefficient U, known also as 



the interaction strength. We note that Eq. ([I]) can be written in a dimensionless form. 
However, since we are interested in the sign effect of the nonlinear coefficient U we keep 
both parameters of Eq. (CQ) in our analysis. 

Equation (CQ) is a simple model which reflects generic features of one-dimensional BEC 
arrays with homogeneous scattering length. The effects of dissipation and inhomogeneities 



are neglected 
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211 ] . Expressions for J and U are known and depend, among 



others, on the lattice spacing, the depth of the optical 
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attice, the s-wave scattering length, 
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2l| . These parameters can be 



and the mass and number of the atoms l[ [13|, 
manipulated in experiments to obtain different desirable configurations of the system. In 
particular the interaction strength U can be tuned through the s-wave scattering length, 
from positive values (repulsive interactions) to negative values (attractive interactions) by 
using either magnetic or laser fields. 



B. Soliton solutions 



In the following we outline the SCQCA for the DNLSE in the spirit of Ref. [22, [28|, |29j, |3( 
to obtain approximate analytical soliton solutions. As mentioned above, the DNLSE is a 
nonintegrable system, i.e. not exact soliton solutions exist. Of course, the solitary solutions 
derived below are only approximations, since SCQCA is only an approximate integration of 



the DNLSE. We note, to the best of our knowledge, that this SCQCA [27|, |28|, |29|, |30[ has not 
been applied to the DNLSE. The importance of this method resides in the fact that it can 
be used without resorting to the knowledge of other nonlinear partial differential equations. 
Therefore it is self-contained. 

In order to proceed with the SCQCA we consider a travelling wave ansatz for an envelope 
complex function reading as 

oo 

tp n (t) = ^2xm(z)exp(im9), (2) 

m=l 

where z = n — voVk t and 9 = kn — e^Ek t+8. Here, k is the quasimomentum, Vk is a velocity, 
Ek is the particle energy, 5 is a phase, and both vq and eo are constants. 

By using the full Taylor expansion of the function (ip n ±i — > exp(±d n )ip n ) , Eq. (TTJ is 
transformed into an operator equation, 

td t Mt) + 2J cos(d n )^ n (t) - W n (t) = 0, (3) 

where W n (t) = U\ip n (t)\ 2 ijj n (t). Inserting the ansatz Eq. (jSj) in Eq. and Fourier trans- 
forming the resultant equation we get 

W m (q) = a m (q)x m (q), (4) 

where 

a m (q) = 2Jcos(km + q) + me Ek + qv v k (5) 
Here tilde marks the Fourier transformed function, e.g. 

oo 

W m (q) = -L / dze* z W m {z), (6) 



where W m is an abbreviation which collect all products of envelope functions that belong to 
the same harmonic e imd , i.e. 

W m (zy me = W n [J2 Xi{z)eA . (7) 

m \ I J 

The SCQCA consist in a formal solution of Eq. (Tj0) for W m (q) and in an expansion of 
the fraction a m (q) (Eq. (jSJ)) for small q [27|, l28(, i.e. a m (q) ~ ^2 n a nm q n . If we consider the 
expansion of the fraction a m (q) up to second order in q, forcing the first order term of the 



expansion to be zero by setting vq = 1, and transforming back to the position space, we find 
a second-order differential equation for the first harmonic \ii namely 

O01X1 - «2i^xi - 3U\xi\ 2 Xi = 0. (8) 
Integrating by parts Eq. (jSJ) and integrating again leads to a solution for the envelope: 



■.B, x _ o. Jcos(k)(e - 1) , / rx rr 7T \ M 



= 2 \l 3^7 " sech ( ^ 2 V(eo - l)z) e lt> , (9) 



for JU cos(k) < and eo > 1 (bright soliton), 



tf(z) = ^ f^fj 1} tanh (VT^Toz) e», (10) 

for JU cos(k) > and eo < 1 (dark soliton). From the SCQCA we obtain also that 
Ej. = — 2J cos(k) and = 2Jsin(fc). Notice that the range of high soliton velocity (rapid 
solitons) corresponds to the values \k\ ~ tt/2. 

Here it is important to remark that Eqs. (Q and/or fflQj) satisfy the so called staggering 
transformation tp n — > ( — l) n ipn (* stands for complex conjugation) [lf|, which has been 
widely used in the literature to consider the cases where the nonlinear coefficient U changes 
the sign. This transformation corresponds in Eqs. ([9]) and/or ( TlOl) to the case k — * 7r — k. 
The solutions that emerge after the transformation are similar to Eqs. (Q and/or (flQj) 
but with the sign in front of U changed. The conditions for the bright and dark regions 
are also interchanged, i.e. JUcos(k) — > —JU cos(k). Notice, however, that with this 
transformation it is not possible to change from the bright region to the dark region or vice 
versa [equivalently, it is not possible to transform Eq. ([9]) into Eq. (flOl) or vice versa]. So, 
in the present case, it is not appropriate to infer conclusions from one region based on the 
information of the other region by simply using this transformation. 



It is worth mentioning that the bright soliton solution, Eq. ([9 



is similar to that fol- 



1 Of ] only when cos(fc) = 1 



lowing from the perturbed Ablowitz-Ladik equation 0, fl B a 
in the amplitude of Eq. (Q. On the other hand, a dependence of the amplitude on the 
quasimomentum k was obtained for bright solitons in Ref. but only for the case bright 
case (JU cos(fc) < 0) in Eq. (JTJ. In the case of dark solitons, this dependence was obtained 
in Ref. [ijj]. So far, the dependence of the soliton amplitudes on the function cos(fc) has not 
been analyzed in relation with their stability. 
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FIG. 1: Normalized probability-density maximum , p(t)/p(0) \p(t) = max(|V^(t)| 2 ]), of the bright 
soliton vs. the time t: (i) k = (or 7r), (ii) k = 2tt/5 (or 3n/5), and (iii) k = 0.497T (or 0.5l7r) with 
J = 1 and \U\ = 1. 

It is also good to keep in mind that a single dark soliton, Eq. (fTUl) . corresponds to the 
case where the sites of the BEC array are equally filled except in the region where the hole 
excitation is moving. However, one can construct pulse solitary waves in the dark region by 
subtracting two identical dark solitons separated some distance d > (see Figs. [5] and [6]), 
i.e. 

^n = ^ +d/2 (t)-^. d/2 (t). (11) 

In the following we shall call u n the dark pulse. Notice that the width of u n depends not 
only on k but also on d. By choosing suitable values for d and eo the form of dark pulses u n 
can become very similar (but not equal) to the form of bright solitons. 



III. INSTABILITY EFFECT 



In this section we first present a conjecture about how to estimate the instability effect 
on the moving solitons. Afterwards, with the help of numerical simulations we compare 
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24| and 



predictions from the well-known standard modulation stability analysis [?], |22|, 
from our conjecture. Special emphasis is done in the analysis of the dark region where 
the standard modulation stability analysis fails to predict analytically the instability of the 
system. We note that the instability of staggered dark solitons was numerically studied in 



Ref. 



16|- 



Here we conjecture that the amplitudes of the approximate soliton solutions, Eqs. ([9]) 
and (flOl) . contain qualitative information of the MI strength regardless of the soliton form. 
Since eo and k are the free parameters of the system, they govern the MI strength on the 
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FIG. 2: Probability density of an infinity train of bright solitons (solid line), i.e. p = | ■ ipi^ (z — 
jL)\ 2 , compared with a modulation plane- wave sin 2 (<5) (dashed line). Here Q = W (eo — 1) and 
eo = 1.1. The amplitudes are normalized to one. 

solitons. So, from a simple inspection of the amplitudes in Eqs. (jSj) and ffTUl) we observe that 
the simplest common factor in both solutions containing eo and k regardless of the signs is 

V = v/|cos(fc)(e -l)|. (12) 

The absolute value in Eq. (112p is used for convenience. Since r\ is proportional to the soliton 
amplitudes, we suppose that the MI strength regardless of the region [bright (J U cos(A;) > 0) 
or dark {JU cos(fc) < 0)] is also proportional to 77. 

Notice that the magnitude of 77 reduces as \k\ — > 71/ 2 vanishing exactly at \k\ = 71/ 2 
(maximum soliton velocity). So, here we can conjecture that MI for moving solitons is present 
in the whole first Brillouin zone except perhaps for \k\ = tt/2, i.e. \k\ G [0, 7r/2) U (tt/2,tt]. 
Moreover, since rj increases as | cos(k) \ — > 1, we can conjecture also that the MI strength for 
moving solitons may increase and reach its maximum at \k\ = 0,tt (breather solutions) and 
may become weak as |eo — 1| — ^ (wide solutions). On the other hand, the MI effect on 
moving solitons can be expected regardless of the JU cos(fc) sign, since both approximate 
bright and dark soliton solutions have been derived. 

The effect of MI on the small-amplitude soliton shapes [Eqs. (jSJ) and (IT0]) ] is of "self- 
defocusing" nature, i.e. a broadening of the soliton width accompanied by exponential-like 
decay of the amplitude can be observed. As an example, in Fig. [TJ it is shown the evolution 
of the normalized probability-density maximum of bright solitons for different values of the 
quasimomentum k. This figure shows that the MI strength can be characterized by studying 
the growth rate of the soliton amplitudes. Since in the present theory only defocusing MI 
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FIG. 3: (Color online) Growth rate for bright solitons, [Eq-©]- a: Versus £q — 1 for (i) k = 

0(or7r), (ii) k = 37r/10 (or 77r/10), (in) k = 2tt/5 (or 3ir/5). b: Versus k for (i) eo = 1.01, (ii) eo = 

1.05, (Hi) €q = 1.1. Numerical estimation T(T/2) with T = 10 (solid line) and analytical estimation 

ki9 ! (£'q) (dashed line). k\ = 0.4 to fit numerical results. 

0.035 

0.030 

iSO.025 
a 

^0.020 
|0.015 
dO.010 
0.005 



/l 


0.030 




// 

®/, 


S0.025 

ca 




// 


^.020 

JS 




It A 


"^0.015 
o 






O0.010 






0.005 






0.000 





(a) 



0.000 

0.00 0.05 0.10 

€ -l 



(b) 



n/2 ii 
k 



FIG. 4: (Color online) Numerical estimation of the Growth rate T(T/2) with T = 10 (solid line) in 
the dark region for u test = it/j^ [i is the imaginary unit and ip^ is given in Eq.©]. The same cases 
as in Fig. [3] are considered. In the panel (a) the analytical estimation for bright solitons K2Q(Eq) 
(dashed lines) is also plotted with K2 = \[2k\ (see caption Fig. [3]). 



effect is observed, the growth rates for moving solitons are decaying 



We note that solitary waves with higher amplitudes than those of Eqs. ([9]) and ( jl~0l) can 
undergo other MI effects as, e.g., self-trapping effect [9j,Ll3| and strong oscillatory instabilities 

In order to check our conjectures above we proceed to estimate a growth rate from the 
standard modulation stability analysis of plane waves in one-dimensional arrays 
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24 



23|, 



251 ] . Here, one looks at the dispersion relation of a phonon wave whose amplitude and 



phase are perturbed by a modulation plane-wave. This yields a dispersion relation that for 
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decaying small- amplitude waves read as [22], |23j, [25] 



tt Q = 2Jsin(fc)sin(Q) - V8J x 



V^l cos(fc) sin 2 (Q/2) + 2Jcos 2 (k) sin 4 (Q/2), 



(13) 



where Qq and Q are the frequency and wavenumber, respectively, of the modulation plane- 
wave, ipo is the amplitude of the phonon wave. The growth rate T of the modulation wave 
can be estimated from the imaginary part 3(f2Q) of Eq. (fT3"|) . i.e. T ~ Q^fig). Notice that 
3(f2 Q ) ^ if JUcos{k) < (bright region) and Q G (0, 2 arcsin[ v /-[/?/'o sec(fc)/(2J)]). 
At first glance, the comparison between a phonon wave modulated by a plane-wave and 

and fflUj) may appear "strange". However, both 



the single soliton solutions in 
problems are strongly linked |7|, |22|, |23j, |24j • Notice that the modulated phonon wave in 
Eq. ( {TBI) can be seen as infinity train ^ = ^2jipn( z ~ jL) of bright solitons, where L is a 
distance between the solitons. In Fig. [2] an example of this comparison is shown. Similar 
construction can be done in the dark region with the help of an infinity train of dark pulses 
[see Eq. (fTTj) ] to emulate a modulated phonon wave. 

A simple estimation of Q(Qq) for bright solitons [Eq. ([9])] can be done by assuming in Eq. 



( TT3l) that ipo is the soliton amplitude [ipo = cos(k)(eo — l)/(— 3£/)] and Q is proportional 
to the inverse width [Q = sj (eo — 1)]. We note that the approximate estimation of Q can 
be done by comparing the envelope shape of a half-cycle oscillation of the modulation wave 
with the shape of the bright-soliton envelope (see Fig. [2]). 

Notice that the case 3(fig) = (dark region: U cos(fc) > 0) means that plane waves 
are stable. However, if one tries to extend this conclusion for solitons, one would wrongly 
conclude that dark solitons are stable. This wrong notion would imply that MI effect is 
negligible in the dark region, which contradicts the fact that approximate dark solitons 
can be systematically derived. Moreover, instability for staggered dark solitons have been 



numerically predicted t 
analysis, Eq. (Tl3l) ^ 5 



16| . Hence, we can conclude that the usual modulation stability 



24j . is not appropriate for deriving conclusions on the stability 
of solitons in the dark region. 

It has been argued that the staggering transformation (k — > n — k) in Eq. (TTBl may 
overcome this problem. However, as mentioned above, by using this transformation the 
conditions for the bright and dark regions are also automatically interchanged. In other 
words, with this transformation one changes from one bright region to other bright region. 



10 



So, in the present case the staggering transformation in combination with the usual mod- 
ulation stability analysis cannot analytically describe MI effect on moving solitons in the 
dark region. 

The remaining question is: how the MI effect for moving solitons can be analytically 
described in the dark region? In order to investigate this question, in the following, we 
numerically estimate the growth rate associate with the MI strength for both the bright 
and the dark regions. Afterwards, we compare numerical results with estimations from the 
standard modulation stability analysis and from the the term 77 in Eq. (TP2|) . We note that 
the analysis in the bright region is presented only for completeness, and in order to compare 
with results obtained in the dark region. 

From numerical simulations in Fig. [T] we can observe that the growth rate depends on 
the time, r = T(t), in a non-trivial form. In fact, the exponential-like nature of the MI effect 
observed in the probability density, p n = |^| 2 , can be described by a exponential function 
of the form p n ~ exp(— 2 f Vdt'). 

The growth rate can be estimated in simulations by using the expression T(t) = 
-dp m ax(t)/(2p max (t)dt), where p max {t) = max(p n (t))dt' /T is an average in time of 

the maximum of p n . This average is used to smooth out small oscillations of p n due to 
the discreteness of the system and is valid for t > T/2. Here T is a small time scale, 
i.e. 2/T ^> r(T/2). Notice that these small oscillations in time of the solitons have been 
discussed in Ref. jll]. 

Since the analytical estimation of the growth rate [r ~ S^q)] can be done only for 
the initial soliton form, we perform a comparison only for a short initial time scale, namely 
t = T/2. 

In Fig. [3] we show the numerical estimations of the initial T value, i.e. T(T/2), for the 
bright region . A comparison with the analytical from the standard modulation analysis 
is also presented. We observe that except for a constant factor (/ti = 0.4) the analytical 
approximation ki^s(Eq) describes qualitatively well the behavior of the growth rate. In 
particular, in Fig. [3k it is possible to observe that the strength of the MI effect reduces as the 
soliton width (amplitude) increases (decreases), i.e. eo — 1 — > 0. In Fig. [3b estimations of the 
growth rate for the first Brillouin zone are presented, showing that MI effects become weak 
as I cos(k)\ — > 0, but they do not completely vanish as predicted by the theory. However, 
T tends to be negligible for very small amplitude solitons, i.e eo — 1, k ~ tt/2. Notice that 
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qualitative predictions from 77 [Eq. (fT2l) ] are identical to the standard modulation analysis 
in the present case. 

Now we consider the case of MI effect on solitons in the dark region (JU cos(k) > 0), 
where Eq. f[T3"l) does not provide any information. However, since we have found approximate 
dark soliton solutions, we should expect also MI effects here. Besides, since 77 [Eq. (TT21) ] is 
equal in both regions, we can expect that the behaviour MI strength on pulse solitons in the 
dark region is similar to that in the bright region. 

Notice that in order to obtain more information about the MI effect on solitons, it is 
convenient to compare both the bright and dark regions. An exact quantitative comparison of 
a bright soliton [Eq. Q] with a dark pulse, u n [Eq. (ITT]) ], is not possible, since their envelopes 
cannot not be exactly equal. So, for the sake of comparison and before we characterize the 
dark pulses u n , we define a sech pulse in the dark region as u^ st = itp^ . Here i is the 
imaginary unit and ip^ is given in Eq.Q. This pulse is identical to a bright soliton, except 
that the relation J U cos(k) > can be considered. 

In Fig. Hh and 0b we show numerical estimations of T vs. eo — 1 and k, respectively, 
for u^ st . The same cases as in Fig. [3] are considered here too. We observe that Fig. 0] 
is qualitatively similar to Fig. [3j However, the T values in the dark region (Fig. 0|) have 
a factor \/2 larger than in the bright region (Fig. [3]). At first glance, it is a surprising 
result, taking into account that plane waves in the dark region are stable. However, it is 
straightforward to observe that 77 ~ Ab/2 or 77 ~ Ad/\/2, where Ab and Ac are the bright 
[Eq. ([9])] and dark [Eq. (ITUl) ] amplitudes, respectively. It means that for identical bright 
and dark pulses with equal amplitude {Ab = Ad), as in Figs. [3] and HI the MI strength, 
represented by 77 is a/2 times higher in the dark region than in the bright one. It also means 
that in order to obtain in the dark region a quantitative similar MI result than in Fig. [3] 
(bright region), one can define a sech pulse in the dark region but with the amplitude and 
width of the dark soliton [Eq. ffTOl) ]. regardless of the soliton form. This assumption can be 
straightforwardly corroborated by numerical simulations. 

The results above show that the strength of the MI effect in the bright and dark regions 
is different. And this difference cannot be estimated with the usual standard modulation 
stability analysis in combination with the staggering transformation. Moreover, it shows that 
the approximate soliton solutions, in fact, do contain information regarding their modulation 
instability. To the best of our knowledge, this type of analysis for moving solitons in the 
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FIG. 5: (Color online) Snapshots of the probability density p n of a Bright (left panel) and dark 
(right panel) breathers at different times [t = (solid line), 200 (dashed line), and 400 (dotted 
line)]. J = 1, U = —1, |eo — 1| = 0.01, k = (left panel), and k = ir (right panel). 
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FIG. 6: Snapshots of the probability density p n of a collision between a bright soliton [Eq. ([9])] and 
a dark pulse [see Eqs. (fTTj) and (fl~0|) ] at different times (top-bottom: t = 0,200, and 400). Arrows 
indicate direction of motion: "— >" for the dark pulse and "<— " for the bright soliton. k = 0.49-/T 
(bright), k = 0.517T (dark), e = 1.02 (bright), e = 0.99 (dark), J = 1, and U = -1. 

dark region and in the whole first Brillouin zone has not been reported, so far. 

For completeness in Figs. [5]and[6]snapshots of the probability-density evolution of solitons 
[Eq. are dark pulses [Eq. (ITT]) ] are shown. In particular, in Fig two breather cases 
(| cos(/c)| = 1) for the bright (Fig. [5^) and dark region (Fig. \5jp) are shown, we observe the 
usual broadening of the pulse and the decaying of the amplitude at different time scales. 
It is important to remark that in the case of the dark region (Fig. [5}d) the collapse of the 
breather starts from the lateral sides, and then moves to the center. Notice that the breather 
maximum does not decay immediately but when the collapse process reaches the center. 

In Fig. [6] an example of a collision between a high- velocity bright soliton and dark pulse 
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is shown. We observe that those high- velocity solitons are stable against collisions and they 
are less prone to MI effect , as expected from the analysis done for [31 and HI 

For very long time scales solitons in Fig. undergo similar collapse as for the breather 
case in Fig. We note that the small amplitude radiation observed in Fig. [H] have been 
investigated with great detail in Refs. [5|. 

Since, the solitons from Eqs. © and (flQj) undergo only "self-defocusing" instability, we 
can conclude that these solutions pose an analytical upper boundary for "defocusing-like" 
solitary waves of the DNLSE. So, solitary waves wit 
instabilities as the oscillatory or the self-trapped [9L Il3l. Il7|. 

Finally, we note that the soliton solutions in Eqs. (Q and fllOp vanish at \k\ = tt/2. This 
value, known as the "zero-dispersion" point, has been examined with great detail in Ref. 



i higher amplitudes can undergo other 



IV. CONCLUSIONS 



Motivated by the problem of coherent matter wave transport in BEC arrays, within 
the first-band approximation, we have studied anew the problem of modulation instability 
(MI) of small-amplitude solitary waves in the discrete Schrodinger equation (DNLSE). For 
that we have developed a self-contained quasicontinuum approximation (SCQCA) for the 
DNLSE to derive approximate analytical soliton solutions. We have used the well-known 
notion that solitons can be considered as a signature of the MI to conjecture that analytical 
soliton solutions following from the approximate integration of the system contain already 
qualitative information of the MI strength on the solitons when propagating. We have 
shown with the help of numerical simulations that this conjecture describes qualitatively 
well the "self-defocusing" MI effect of small-amplitude solitons in the dark region where the 
standard modulation stability analysis of planewaves and/or the staggering transformation 
do not provide any information. Though planewaves are stable in the dark region, we have 
shown that for identical pulses in the bright and dark region the strength of the "self- 
defocusing" MI is higher in the dark region than in the bright one. This fact was shown not 
only numerically but also by following the conjecture posed above. 

Since the soliton solutions derived here only present "self-defocusing" instability, their 
amplitudes can be considered as analytical upper boundaries for this instability. 



14 



Last but not the least, the analysis proposed here can be straightforwardly extended to 
higher spatial dimensions (work in progress) where the motion of solitary waves and vortices 
can be also observed and analyzed. 
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